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Abstract 

We investigate the uncertainty in the strong coupling as(M§) when allowing it to be a free 
parameter in the recent MMHT global analyses of deep-inelastic and related hard scat¬ 
tering data that was undertaken to determine the parton distribution functions (PDFs) 
of the proton. The analysis uses the standard framework of leading twist fixed-order 
collinear factorisation in the MS scheme. We study the constraints on ag(M|) coming 
from individual data sets by repeating the NNLO and NLO fits spanning the range 0.108 
to 0.128 in units of 0.001, making all PDFs sets available. The inclusion of the cross sec¬ 
tion for inclusive tt production allows us to explore the correlation between the mass nit of 
the top quark and ag(M|). We find that the best fit values are a${M (|) = 0.1201 ±0.0015 
and 0.1172 ± 0.0013 at NLO and NNLO respectively, with the central values changing to 
as(Af|) = 0.1195 and 0.1178 when the world average of ag(M§) is used as a data point. 
We investigate the interplay between the uncertainties on as(M^) and on the PDFs. In 
particular we calculate the cross sections for key processes at the LHC and show how the 
uncertainties from the PDFs and from as(M^) can be provided independently and be 
combined. 


1 Introduction 

There has been a continual improvement in the precision and in the variety of the data for 
deep-inelastic and related hard-scattering processes. Noteworthy additions in the years since 
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the MSTW2008 analysis |T] have been the HERA combined HI and ZEUS data on the total [2] 
and charm structure functions [3j, and the variety of new data sets obtained at the LHC, as 
well as updated Tevatron data (for full references see |4]). Moreover the procedures used in 
the global PDF analyses of these data have been refined, allowing the partonic structure of the 
proton to be determined with ever-increasing accuracy and reliability. These improvements 
are important as it is necessary to quantify the Standard Model background as accurately as 
possible in order to isolate possible experimental signals of New Physics. One area that needs 
careful attention, at the present level of accuracy, is the treatment of the strong coupling, as 
itself, in the global analyses. Here we extend the recent MMHT2014 global PDF analysis [3] to 
study the uncertainties on as and their implications for predictions for processes at the LHC. 


2 Treatment of in the MMHT2014 analysis 

We refer to Fig. [T| for an overview of the treatment and of the values of as obtained in the 
MMHT2014 global PDF analysis [4|. At both NLO and NNLO the value of as(M|) is allowed 
to vary just as another free parameter in the global fit. The best fit values are found to be 

a s,NLo(^l) = 0.1201 (1) 

«s,NNLo(Aff) = 0.1172, (2) 

as indicated by the dark arrows in Fig. [l] The corresponding total y 2 profiles versus as(M§) 
are shown in Fig. [2j These plots clearly show the reduction in the optimum value of as(M§) 
as we go from the NLO to the NNLO analysis. In the next section we show how the individual 
data sets contribute to make up this y 2 profile versus as(M|). 

It is sometimes debated whether one should attempt to extract the value of cts(Mf) from 
the PDF global fits or to simply use a fixed value taken from elsewhere - for example, to use the 
world average value [5J. However, we believe that very useful information on the coupling can 
be obtained from PDF fits, and hence have performed fits where this is left as a free parameter. 
As the extracted values of as(M§) in the NLO and NNLO MMHT2014 analyses [Tj reasonably 
bridge the world average of a s (M|) = 0.1185 ± 0.0006 |>], we regard these as our best fits. We 
note it is a common result in PDF analyses, and in other extractions of the strong coupling, for 
the best fit value to fall slightly as the order of the theoretical calculations increases. However, 
in order to explore further, as well as leaving as(M^) as a completely independent parameter, 
the MMHT2014 analyses were repeated including the world average value (without the inclusion 
of DIS data to avoid double counting) of as(M^) = 0.1187 ± 0.0007 as a data point in the fit. 
This changed the preferred values to 

«s,NLo(d7f) = 0.1195 and as,NNLo(Aff) = 0.1178, (3) 

as indicated by the short arrows in Fig [TJ Each of these is about one standard deviation away 
from the world average, so our PDF fit is entirely consistent with the independent determina¬ 
tions of the coupling. Moreover, the quality of the fit to the data (other than the single ‘data’ 
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Figure 1: The dark arrows indicate the optimal values of as(Mf) found in NLO and NNLO fits 
of the MMHT2014 analysis [4j. The dashed arrows indicate the values found in the MSTW2008 
analysis [T], We also show the world average value, which we note was obtained assuming, for 
simplicity, that the NLO and NNLO values are the same - which, in principle, is not the case. The 
short arrows are also of interest as they indicate the NLO and NNLO values which would have 
been obtained from the MMHT2014 global analyses if the world average value (obtained without 
including DIS data) were to be included in the fit. However, the default values o^nlo = 0.120 and 
«s,nnlo = 0.118 were selected for the final MMHT2014 PDF sets ‘for ease of use’; indeed, the 
small values of Ay 2 are the minute changes in y 2 lobal ' n S°' n g from the optimal to these default 
fits. 


point on as(Mf)) increases by about 1.5 units in y 2 at NLO and just over one unit at NNLO 
when the coupling was included as a data point. 

However, ultimately for the use of PDF sets by external users it is preferable to present 
the sets at common (and hence ‘rounded’) values of «s(M|) in order to compare and combine 
with PDF sets from other groups, for example as in [6, 0 8) [9]. At NLO in the MMHT2014 
analysis [4] we hence chose as(M z ) = 0.120 as the default value for which the PDF sets with 
full error eigenvectors are made available. This is essentially identical to the value for the best 
PDF fit when the coupling is free, and still very similar when the world average was included 
as a constraint. At NNLO, as(Mf) = 0.118 was chosen as a rounded value, very near to 
both the best fit value and the world average, and the fit quality is still only 1.3 units in y 2 
higher than that when the coupling was free. This is extremely close to the value determined 
when the world average is included as a data point. Hence, in MMHT2014 [4], we chose to 
use ag(Mf) = 0.118 as the default for NNLO PDFs, a value which is very consistent with the 
world average. At NLO we also made a set available with ag(Mf) = 0.118, but in this case 
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Global fit: Xo = 3268.6 for 3040 pts. 



as(Afi) 

Global fit: Xo = 2719 for 2707 pts. 



Figure 2: The upper and lower plots show total x 2 as a function of the value of the parameter 
as(M'z) for the NLO and NNLO MMHT2014 fits, respectively. 


the x 1 2 increases by 17.5 units from the best fit value. In [Tj we also made available PDF sets 
corresponding to the best fit for cts(Mf) values ±0.001 relative to the default values in order 
for users to determine the as (Mf)-uncertainty in predictions if so desired. We will return to 
the issue of PDF ±cts(Mf) uncertainty later. 

Before we continue we should specify how the running of as(Q 2 ) is treated. There is 
more than one definition of the coupling commonly used in QCD phenomenology. Although 
the various prescriptions are all formally equivalent since they differ only at higher orders, 
numerical differences of the order of up to 1% can occur. We use the definition based on the 
full solution of the renormalisation group equation, in MS scheme, at the appropriate order, 
with boundary condition defined by the value of as(M^). This is identical to the definition 
in public codes such as PEGASUS [10] and hoppet na, and is now effectively the standard in 
PDF analyses]]] It differs, for example, from solutions to the renormalisation group equations 
truncated at a particular order. 

1 In MS scheme this involves discontinuities at flavour transition points at NNLO. For a suggestion for a 

smooth transition in a physical scheme see m- 
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3 Description of data sets as a function of as 


The NNLO MMHT2014 global analysis [1] was based on a fit to 40 different sets of data 
on deep-inelastic and related hard scattering processes. There were 10 different data sets 
of structure functions from the fixed-target charged lepton-nucleon experiments of the SLAC, 
BCDMS, NMC and E665 collaborations, six different neutrino data sets on F 2 , xF 3 and dimuon 
production from the NuTeV, CHORUS and CCFR collaborations, two Drell-Yan data sets from 
E886/NuSea, six different data sets from HERA involving the combined HI and ZEUS structure 
function data, seven data sets from the Tevatron giving the measurements of inclusive jet, W 
and Z production by the CDF and DO collaborations and, finally, nine data sets from the 
ATLAS, CMS and LHCb collaborations at the LHC. In addition, the NLO fit also used jet 
data from the ATLAS, CMS and HI and ZEUS collaborations, which were not used at NNLO 
because it was judged that at present there is not sufficient knowledge of the full jet NNLO 
cross section; jet production at the Tevatron, on the other hand, is much nearer to threshold 
than at the LHC, so the threshold approximation to the full NNLO calculation is much more 
likely to provide a reliable estimate in this case. The goodness-of-fit quantity, x 2 ) for each of 
the data sets, n = 1, ...40, is given for the NLO and NNLO global fits in the Table 5 of |4], and 
the x 2 definition is explained in Section 2.5 of the same article. The references to all the data 
that are fitted are also given in [4], 

In the NNLO global fit of [4], let us denote the contribution to the total x 2 from data set 
n by Xn O' Here we explore the x 2 profiles as a function of «s(M|) by repeating the global 
fit for different fixed values of «s(Mf) in the neighbourhood of the optimum value given in 
(j2j. The results are shown in Figs. 3-7, where we plot the x 2 profiles for each data set n as 
the difference from the value at the global minimum, x 2 o> when varying as(Mf). The points 
(•) in Figs. 3-7 are generated for fixed values of as(Mf) between 0.108 and 0.128 in steps of 
0.001. These points are then fitted to a quadratic function of as(M|) shown by the continuous 
curves. By definition we expect the profiles to satisfy (x 2 — Xn o) — 0 at as(M^) = 0.1172, 
corresponding to the value of as(M^) at the NNLO global minimum. Ideally, a data set should 
show a quadratic minimum about this point. Of course, in practice, the various data sets may 
pull, in varying degrees, to smaller or larger values of as(Mf). There is a small amount of 
point-to-point fluctuation for the values of (y 2 ~ Xn o)> even near the minimum, but near the 
minimum this is generally only at the level of fractions of a unit in x 2 for a given data set. The 
fluctuations become larger as we go to values of as(M§) far from the minimum, particularly 
for lower ots{M §), mainly because changes in x 2 with small changes in as (Mg) are becoming 
much greater. However, some of the “jumps” for individual sets near cts(Mf) = 0.108 imply 
that the global minimum in y 2 f° r this choice of as(Mf) is rather flat in certain parameter 
directions, with some relatively easy trade-off between the data sets which are poorly fit, and 
a transition to a different, approximately degenerate global minimum occurring with a small 
change in as(M|). Indeed, we have verified that at as(M§) = 0.108 there is a local minimum 
where “jumps” are eliminated, but with slightly higher global x 2 than the result where there 
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are “jumps”. This highlights the fact that the PDF uncertainty is difficult to define properly 
for a value of «s(M§) which is far from optimal and leads to many data sets being badly fit. 

We repeat this exercise at NLO. Then the profiles will satisfy (\n ~ Xn o) = 0 at as(M^) = 
0.1201. We include in the plots the NLO points (as triangles) and show the corresponding 
quadratic fit by a dashed curve. In Fig. 8 we show the \n profiles for the LHC and HERA jet 
data that were included in the NLO fit. Here the bullet points and profile curve correspond to 
the NLO fit. These data were not included in the NNLO fit. 

The fixed-target structure function data in the first 14 plots in Figs. 3 and 4 have been 
available for several years. These data play an important role in constraining the value of 
There is some tension between these data sets. The BCDMS (and also the E665) data 
prefer values of «s(M|) around 0.110. On the other hand the NMC data prefer values around 
0.122; and the SLAC F%' d data prefer as(M f) values around 0.115 and 0.122 respectively. The 
neutrino F 2 and xF 3 data prefer as(M§) ~ 0.120; while neutrino dimuon production has little 
dependence on «s(Mf), since the extra B(D —>■ p) branching ratio parameter (see eq.(19) of 
m can partially compensate for the changes in as(M §). 

The NNLO corrections to the structure functions are positive and speed up the evolution, 
leading to smaller optimum values of as(Mf) than those at NLO, such that the spread of 
optimum values of ag(M|) for the different data sets is somewhat reduced. Thus the overall 
fit to this subset of the data is marginally better at NNLO. The difference «s,nnlo < «s,nlo is 
clearly evident in the majority of the corresponding plots. 

The recent combined HI and ZEUS structure function data from HERA prefer a value of 
as {Mg) of about 0.120 at NNLO. Perhaps the only surprising result is the as(M§) behaviour 
of the combined data for F^ mrm which prefers a very low value of as(M^) at NNLO, whereas 
the uncombined data had a perfect quadratic behaviour about 0.118, see Fig.5 of ra- Note, 
however, that the combined data contains some points at the lowest Q 2 which were not available 
as an individual data set. These data, particularly at low Q 2 , are sensitive to the value of the 
charm mass m c , and there is a correlation between its value and a s (M|) [Hj. This will be 
studied again with the up-to-date data in a future article. 

The longitudinal structure function F Iy leads off with an as term, and so the value of 
(Xn ~ Xn, o) depends more sensitively on ag(Mf). The NNLO plot shows an excellent quadratic 
dependence on ag(Mf), centred at 0.118. The NNLO coefficient functions for Fl(x, Q 2 ) (15, 16] 
are positive and significant, and the NLO fit tries to mimic these with a higher value of «s(M|). 
Indeed, the data for F L , and also the E866/NuSea pp Drell-Yan cross sections data, are clearly 
more quadratic at NNLO than at NLO, with minima closer to the best-fit values. This indicates 
a strong preference for the NNLO description, which is not so apparent if only the global best- 
fit values xl,o are known. As the E866/NuSea data for pd/pp Drell-Yan production are a ratio 
of cross sections, the sensitivity to the value of as(Mf) is small. 

The Tevatron data, as well as the ATLAS W ± ,Z production data and the ATLAS high- 
mass Drell-Yan data, show, at NNLO, a$(M §) profiles with quadratic behaviour with minima 
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Figure 3: Xn profiles obtained when varying as(M|) for the subset of data from deep-inelastic 
fixed-target experiments. The results from the NNLO global fits are shown by bullet points 
(and a continuous curve), while those from the NLO global fits are shown by triangles (and a 
dashed curve). The plots are continued in next figure. 
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Figure 4: \n profiles obtained when varying as(M^), for the subset of data from deep-inelastic 
fixed-target experiments. The results from the NNLO global fits are shown by bullet points 
(and a continuous curve), while those from the NLO global fits are shown by triangles (and a 
dashed curve). (Continued from the previous figure.) 
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Figure 5: Xn profiles obtained when varying cts{M f) coming from the Drell-Yan fixed-target 
experiments and from the combined HI and ZEUS measurements at HERA. The results from 
the NNLO global fits are shown by bullet points (and a continuous curve), while those from 
the NLO global fits are shown by triangles (and a dashed curve). 
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Figure 6: xl, profiles obtained when varying as(Mf), coming from the subset of data of the CDF 
and DO Tevatron experiments, together with the plot for the ATLAS W and Z production data. 
The results from the NNLO global fits are shown by bullet points (and a continuous curve), 
while those from the NLO global fits are shown by triangles (and a dashed curve). 
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Figure 7: xl, profiles obtained when varying ag(M|), from the subset of data collected by the 
LHC experiments. The results from the NNLO global fits are shown by bullet points (and a 
continuous curve), while those from the NLO global fits are shown by triangles (and a dashed 
curve). The Xn profiles for tt data are shown in Fig. [ 9 ] and discussed in Section [ 4 J 
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Figure 8: Xn profiles for jet data sets, included in the NLO fit, but not in the NNLO fit, when 
varying a s (M§). 
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close to the best fit values. Again, the profiles are improved to those at NLO. The counter 
example are the LHCb data, which have profiles which are more reasonable at NLO than at 
NNLO. In general, the charge-lepton asymmetry measurements arising from W ± production 
at the Tevatron and the LHC, which are a ratio of cross sections, have much less constraint on 
the value of a,s(M§). 

Judging from the values of (y 2 ~ Xn o) aw ay from the different minima of the various data 
sets or, rather, the steepness of the quadratic forms in a s (M|), we see that there is a tendency 
for data at lower energies or lower Q 2 to have more constraint on the optimum global value of 
as(M'z). This is to be anticipated, as we will see in Section |b} 


4 tt data: m t — as correlation 

There is a particularly strong, but also complicated, relationship between the value of as(M^) 
and the fit to data on the inclusive cross section for tt production, a t t- We show the y 2 profiles 
at NLO and NNLO in Fig. [ 9 } Clearly there is a preference for a lower value of as(M^) at NNLO 
than at NLO, and a strong constraint in both cases, with y 2 increasing by a large number of 
units, certainly compared to the number of data points, for small changes in as(M f). Indeed, 
nominally a t t provides one of the strongest constraints of any data set for the lower limit 
of «s(Mf) at NLO and the upper limit of as(Mf) at NNLO. However, the picture is more 
complicated than for other data sets due to the very strong correlation with the value of the 
mass m t of the top quark. 

In the global fits the theory calculation of (J f j is performed with a preferred value of the top 
quark pole mass of 172.5 GeV, since this is the default in PYTHIA, used to extract the cross 
section in many of the measurements. Moreover, the majority of the cross sections are quoted 
for this value of m t . This value is also consistent with the world average of the measured value 
of 173.34 GeV [5j. However, we allow a 1 GeV uncertainty on the value of m t which can be 
thought of as accounting for the uncertainty in the value of m t itself and also for the small 
variation in the extracted cross sections with rri t used; in general this is about a third the size 
of the variation of the calculation of a t t with mt, and the net effect is an effective uncertainty a 
little lower than 1 GeV. To be specific, m t is left as a free parameter in the fit, but there is a y 2 
penalty of y 2 = (m t — 172.5 GeV) 2 applied to keep the value close to the preferred value. This 
penalty is included in the values in Fig. |9j The allowed variation in rn t away from the preferred 
central value of 172.5 GeV results in the NLO fit preferring a low value of mt = 171.7 GeV and 
the NNLO fit preferring a high value of m t = 174.2 GeV. The low value of mt in the global fit 
and the high value of «s(M|) preferred by <j t t when «s(Mf) is varied, both occur for the same 
reason. That is, the NLO cross section tends to undershoot the data, and raising «s(Mf) and 
lowering mt both raise the cross section, leading to better agreement. 

The NNLO correction to the cross section in the pole mass scheme is moderate, but large 
compared to the most precise data, and hence the NNLO cross section tends to be too high. This 
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Figure 9: y 2 profiles for tt data in the NLO (left) and NNLO (right) fits, when varying as(M^). 

leads to the opposite pulls to those at NLO, i.e NNLO prefers «s(Mf) low and m t high. Within 
the global fit we find that the allowed variation with accompanying penalty for deviations from 
mt = 172.5 GeV results in mt values at the best fit values of as(Mf) which are of order 1 — 2cr 
away from either our preferred value or the world average, so have no particular inconsistency, 
but it is useful to examine the interplay between as (Mf) and m t in rather more detail. 


4.1 Effect on X | obal 


oiiaii^co 


First we investigate the quality of the global fit as a function of both a,j(M|) and mt. This 
is shown in Fig. 10, where we plot y 2 lobal versus rn t at several different values of as(M|) (In 
these plots m t is varied with no y 2 penalty for deviations away from the “preferred” value). 
At NLO one can see that regardless of rn t the best global fit is always obtained quite clearly 
for as(Mz) near to 0.120, with the fit quality for as(Mf) = 0.119 or as(M f) = 0.121 each 
being a few units worse at all values of m t . ft is only for mt > 180 GeV that the quality for 
«s(M|) = 0.121 approaches that of 0.120 and the best fit would be for a 5 (M|) ~ 0.1205. 
At this mass the global y 2 is about 10 units above the minimum though. Similarly at NNLO 
as(Mz) = 0.117 gives a lower y 2 lobal for all masses between about 166 GeV and 181 GeV, when 
as(M \) = 0.116 and as(Mf) = 0.118 respectively give the same y 2 lobal values. Hence, even 
completely unreasonable variations of ~ 7 — 10 GeV result in changes of the best fit values of 
a s (M|) of only r^j 0.0005. We do note, however, that without a penalty for m t variation the 
best global fits are at m t = 168 GeV and m t = 180 GeV at NLO and NNLO respectively, so 
some penalty is clearly necessary. Ultimately, the value of «s(Mf) determined by the global fit 
is very insensitive to the value of rn t used, and indeed, to the a t i data, because these correspond 
to relatively few data points. Indeed, if these are left out of the global fit the change in the 
optimum value of as(Mf) is only of order 0.0001 — 2 at NLO and NNLO. However, the interplay 
between as(M^) and rn t is more dramatic for the a t t data alone, as we will now show. 
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NLO: Global x 2 vs. m t for various as 



NNLO: Global x 2 vs. m t for various as 



m t [GeV] 


Figure 10: Global y 2 minima as a function of the top mass mt, for different fixed values of 
as {Mg). There is no y 2 penalty for varying mt- 


4.2 Effect on x% to changes of m t and as(M^) 


The equivalent plots to Fig. 10 are shown in Fig. 11 for the fit quality to the inclusive o t t 
cross section data. Again, there is no penalty applied for rn t variation. At NLO it is clear 
that, except for very low values of mt, the best fit is achieved for higher values of as(M|), he. 
a's(Mf) = 0.121 or for rn t > 172 GeV, as(Mf) = 0.122. Indeed, the best possible fit to the top 
cross section data is for rn t ps 172 GeV and ccs(Mf) = 0.122. However, the improvement in Xti 
compared to < 25 (M§) = 0.120 for this mass is only ~ 2 units - far less than the deterioration 
in the y 2 for the rest of the data when going from as(Mf) = 0.120 to 0.122. Overall the 
minimum y 2 achieved for any as{M%) is quite flat with mt, changing by at most 2 units for 
168 GeV < m t < 178 GeV. However, it is clear that the variation of yT is different for different 
values of as (Mg). As as(M §) decreases there is a preference for a smaller mass, hence if the 
central value of rn t had been chosen higher than 172.5 GeV for example, the best fit to a t t 
would be for a higher value of as(Mf). The constraint on as(M^) in the upper direction 
would be weakened slightly, however this data set does provide a significant constraint in this 
direction. If the penalty had been less severe, e.g. a increase in y 2 f for A m t = 2 GeV rather 
than Amt — 1 GeV, the best value of m t and as (Mg) would not change significantly, as the 
fit quality does not improve for masses of m t < 171.7 for any as(M§), even discounting the 
penalty. However, the y 2 f curves f° r lower values of as (Mg), i.e. 0.119 and 0.118 are falling 
quite steeply as rn t decreases in the vicinity of rn t = 172 GeV, so the increase in Xti with 
decreasing as{M f) seen in Fig. [ 9 ] (left) would be less severe if m t was allowed to choose smaller 
values, and the constraint on the lower values of a s (M|) would be reduced somewhat. Hence, 
at NLO, alternative treatments of rn t would allow a slightly higher best fit as(M^) than the 
default treatment, and a little scope for a relaxation of the lower limit on as (Mg). 


At NNLO it is again clear that higher values of as(M^) prefer higher values of m;. However, 
for as{Mz) = 0.118 or 0.119 the value of m t corresponding to the best fit is mt = 180 GeV 
or more. Again, there is little variation in the best value of y 2 f for 168 GeV< rrq < 178 GeV, 
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Figure 11: y 2 values for inclusive tt cross section data at the global minimum, as a function of 
the top mass mt, for different fixed values of as(M§). There is no y 2 penalty for varying mt■ 


but the best fit is achieved for a s (M f) = 0.115 or 0.116 [^J only becoming a s (M f) = 0.117 at 
nit = 178 GeV. In this case if the penalty for variations in rn t away from the default central 
value were relaxed it would make little difference, as even for as(M §) = 0.115 the best fit is 
for mt ~ 172 GeV. It might allow slightly better fits for a?s(M|) ~ 0.110, but this would 
have no influence on the overall constraint on as(Mf), which is constrained by many data not 
to be much lower than 0.115. A potential increase in m t , either by change of default central 
value, or a relaxation of the penalty, would allow for a potentially a slightly higher value of 
rn t for the best fit, as the minimum possible yT almost completely flat between 172 GeV 
< rn t < 176 GeV. This would be accompanied by a slight increase in ets(Mf). It would 
also allow a little relaxation in the constraint on higher values of a,s(M§). The yT curves for 
as(M’z) = 0.118 and 0.119 are decreasing with increasing rn t in the vicinity of rn t = 174 GeV, 
and a higher allowed value of m t would enable the increase in y 2 with as(M^) in Fig. 11 (right) 


to be less steep. Hence, at NNLO alternative treatments of rn t would allow a slightly higher 
best fit value of a?g(M|) than the default treatment, and a little scope for a relaxation of the 
upper limit on as (M|). 

Hence, the overall conclusion is that some added freedom in m t would lead to potentially 


rather small changes in the minima of the y 2 curves in Fig. 11, but a reduced rate of increase 
of y 2 away from the minima. The implications of this will be discussed in the next section. 


5 Uncertainty on and calculation of PDF+Cfs(M|) 

uncertainty 

First, recall that in the MMHT2014 analysis [1] we determined the uncertainties of the PDFs 
using the Hessian approach with a dynamical tolerance procedure. We obtained PDF ‘error’ 

2 Our constraint on ag(M|) is very consistent with that in mi¬ 
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eigenvector sets, each corresponding to 68% confidence level uncertainty, where the vectors are 
orthogonal to each other and span the PDF parameter space. 


In order to determine the uncertainty on ats(M §) at NLO and NNLO we begin by using the 
same technique as in the MSTW study of Ref. [13]; that is, for the ‘error’ eigenvectors we apply 
the tolerance procedure to determine the uncertainty in each direction away from the value at 
the best fit when one data set goes beyond its 68% confidence level uncertainty. The values at 
which each data set does reach its 68% confidence level uncertainty, plus the value of as(Mf) 
for which each data set has its best fit (within the context of a global fit) are shown at NLO and 
NNLO in Fig. 12 However, unlike Fig. 7 of [13] we do not show all data, as with the increased 


number of sets there are now too many to show clearly on a single figure. Moreover, as seen 
earlier, many data sets have very little dependence, and hence produce very little constraint. 
Hence, we show those where both limits are within the range of as(M f) explicitly studied, i.e. 
0.108 — 0.128 or where one limit is within 0.005 of the best fit value of as(M z ). None of the 
data sets omitted using these criteria have a significant pull on ag(M§). 


The dominant constraint on cts/Mf) in the downwards direction at NLO is from the top 
pair cross section data and, using the dynamical tolerance procedure, gives an uncertainty of 
Aas(Mf) = —0.0014. In the upwards direction it is the BCDMSp data with an uncertainty 
of Aag(Mf) = +0.0012. At NNLO the dominant downward constraint comes from NuTeV 
F 3 (x, Q 2 ) data which gives Aag(Mf) = —0.0012 and in the upwards direction it is the top pair 
cross section data, where the uncertainty is Aag(Mf) = +0.0008. 


There are a number of other data sets which give almost as strong constraints. For instance, 
at NLO in the downwards direction we find that SLAC deuterium data give A as(M' z ) = 
—0.0018 and in the upwards direction HI jets give A a s (M z ) = +0.0019. At NNLO in the 
downwards direction SLAC deuterium data and CDF jet data give Aag(Mf) ~ —0.0014, and 
in the upwards direction, at NNLO, the BCDMSp data give Aag(Mf) = +0.0014. In all cases 
there are other data sets that are not much less constraining than those mentioned explicitly. 
Hence, in no case is it a single data set which is overwhelmingly providing the dominant 
constraint on the upper or lower limit of ag(M|). Similarly, no single data sets would change 
the central value by more than 0.001 if it were to be omitted. 


Two of the four dominant constraints nominally come from 0 +, and at NLO we have 
a s (M 2 z ) = 0.1201^;[]^4 and at NNLO a s (M|) = 0.1172^;“^. However, in the previous 
section we highlighted the interplay between ag(M§) and m t when examining the fit quality 
of the cr t i data. We demonstrated that if some extra flexibility is allowed on the choice of 
central value of rn t and/or on the 1-cr uncertainty that is used, then the constraints are re¬ 
laxed to some degree. Hence, we are reluctant to treat the constraint from the data on a t i 
completely rigorously. In order to see quite how we should deal with the constraints nominally 
due to these data we first check which data sets provide the next tightest constraint. If we 
were simply to ignore the constraints from a t t we would find a change in uncertainty at NLO 
of Aas(Mf) = -0.0012 -+ -0.0017 and at NNLO Aa s (M z ) = +0.0008 -+ +0.0014. These 
are significant, but hardly dramatic changes, and it would be no surprise if some alternative 
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treatment of the default top mass resulted in changes of a similar type. Hence, it might be 
suitable to take these values as a simple alternative, arguing that the constraints from a t j are 
not sufficiently greater than those from other data sets either to ignore the possible effects of 
alternative treatments of the mass rn t or to warrant a completely thorough investigation at this 
stagej^] However, there is the additional feature to note - whichever criterion we use we have 
some, albeit not too dramatic, asymmetry in the ag(M f) uncertainty. There is no strong reason 
to apply this slight asymmetry, as the y 2 profile for the global fit follows the quadratic curve 
very well at both NLO and NNLO, and the degree of asymmetry obtained using the dynamical 
tolerance procedure is arguably within the “uncertainty of the uncertainty”. Hence at NLO 
and NNLO we average the two uncertainties (obtained without the a t j constraint) obtaining 

«s, nlo (Aff) = 0.1201 ±0.0015 (4) 

as,NNLo(M|) = 0.1172 ±0.0013. (5) 

This corresponds to A NLO y 2 lobal = 10.3 and A NNLO y 2 lobal = 7.2. These are the sort of tolerance 
values typical of the majority of PDF eigenvectors. 

Each of these values of ag(Mf) is within la of the world average (without D1S data) of 
0.1187 ± 0.0007, though in opposite directions. As noted earlier, the inclusion of ag(Mf) as 
a data point leads to values of 0.1178 and 0.1195 at NNLO and NLO respectively. These 
are somewhat closer to the world average, and very near to 0.118 at NNLO, but still quite 
close to 0.120 at NLOj^] Hence, we interpret the values in eqs. ([I]) and (J5]) as independent 
measurements of ag(Mf), but acknowledge that at NNLO taking both this determination and 
the world average into account a round value of as(M^) = 0.118 is an appropriate one at which 
to present the PDFs. At NLO we would recommend the use of ag(Mf) = 0.120 as the preferred 
value for the PDFs, but have made eigenvector sets available at ag(Mf) = 0.118. If a value of 
ag(M|) = 0.119 were desired the average of the results at ag(Mf) = 0.118 and 0.120 would 
provide an excellent approximation. 

When considering the uncertainty on the prediction for a physical quantity we should include 
the uncertainty on «g(M|), as well as that on the PDFs. This is particularly important for 
cross sections that at leading order are proportional to a power of the coupling, such as a t t or 
a Higgs; which are proportional to a 2 s . A naive procedure would be to compute the error as 

Act = yj (Actpdf) 2 ± (A a as ) 2 (6) 

3 The constraint from cr t i data does provide the dominant constraint in one direction for eigenvector 15 at 
NNLO. However, very nearly as strong a constraint is provided by other data sets and the eigenvector only 
provides at the very most 40% of the uncertainty on one distribution, the gluon, at any x value, in practice 
at high x. Hence, a slightly increased tolerance for this eigenvector would have minimal impact on any PDF 
uncertainties. 

Waking a weighted average of the values in eqs. Q and ([5| would result in values slightly nearer to the 
world average, reflecting the fact that the dynamical tolerance procedure used to determine the uncertainty 
results in a Ay|; lobal > 1. 
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where A a as is variation of the cross section when as(M%) is allowed to vary over a given range. 
However, it is inconsistent to use different values of as in the partonic hard subprocess cross 
section and in the PDF evolution. Moreover, in a global PDF analysis, there are non-negligible 
correlations between the PDFs and the value of as. 

In the MSTW study |13| of the PDF+as(Mf) uncertainties arising from the MSTW2008 
analysis we advocated using our best fit value of as(Mf) as the central value for PDF pre¬ 
dictions, and then provided additional eigenvector sets at ±0.5cx and ±lcr values of as(Mf). 
The uncertainty was then calculated by taking the envelope of the predictions using all these 
eigenvector sets. This still seems like an appropriate algorithm for use with the dynamical 
tolerance procedure of obtaining uncertainties. However, it can only really be applied if the 
central prediction is obtained using the PDFs defined at the best fit value of as (M|), which is 
no longer the case, and moreover was a rather complicated and time-consuming procedure. 

Since the MSTW study [13] was undertaken it has been shown that, within the Hessian 
approach to PDF uncertainties, the correct PDF +as(Mg) uncertainty on any quantity can be 
obtained by simply taking the PDFs defined at «s(M|) ± A«s(Mf) and treating these two 
PDF sets (and their accompanying value of as(M|)) as an extra pair of eigenvectors [18j. In 
short, the full uncertainty is obtained by adding the uncertainty from this extra eigenvector 
pair in quadrature with the PDF uncertainty. So we are back to the naive form ([6]), but 
now, importantly, with the correlations between the PDFs and as included. This has the 
advantages of both being very simple, but also separating out the «s(M|) uncertainty on a 
quantity explicitly from the purely PDF uncertainty. Strictly speaking, the method only holds 
if the central PDFs are those obtained from the best fit when as (Mg) is left free, and if the 
uncertainty Aas(Mf) on ag(M|), that is used, is the uncertainty obtained from the fit. If we 
use PDFs defined at as(M §) = 0.118 at NNLO we are still very near the best fit, and the error 
induced will be very small. At NLO a larger error will be induced by using the PDFs defined 
at «s(M|) = 0.118 than those at ag(Mf) = 0.120. Any choice of Aas(Mf) of 0.001 — 0.002 
should only induce a small error. Hence, overall we now advocate using this approach with 
NLO PDFs defined at as(M §) = 0.120 and NNLO PDFs defined at as(M |) = 0.118. The 
value of Aas(Mf) is open to the choice of the user to some extent, but it is recommended to 
stay within the range Aois(Mf) that we have found. 

In Section[7]we apply the above procedure to determine the PDF+a,s(Mf) uncertainties on 
the predictions for the cross sections for benchmark processes at the Tevatron and the LHC, 
but first we examine the change in the PDF sets themselves with as(Mf). 


6 Comparison of PDF sets 


It is informative to see the changes in the PDFs obtained in global fits for fixed values of 
as(M'z) relative to those obtained for the central value; we only consider the NNLO case here, 


but note that the NLO PDFs behave in a similar way. These are shown in Figs. 13-15 for 
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Figure 12: The upper and lower plots show the value of as(M §) corresponding to the best fit, 
together with the upper and lower lcr constraints on as(M|) from the more constraining data sets 
at NLO and NNLO respectively. 
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Figure 13: Percentage difference in the NNLO gluon and strange quark PDFs at Q 2 = 10 4 
GeV 2 relative to central (as(M§) = 0.118) set for fits with different values of as, with the 
percentage error bands for the central set also shown. 


Up quark (NNLO), Q 2 = 10 4 GeV 2 


Down quark (NNLO), Q 2 = 10 4 GeV 2 



X X 


Figure 14: Percentage difference in the NNLO up and down quark PDFs at Q 2 = 10 4 GeV 2 
relative to central (as(M§) = 0.118) set for fits with different values of as, with the percentage 
error bands for the central set also shown. 
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Up valence (NNLO), Q 2 = 10 4 GeV 2 


Down valence (NNLO), Q 2 = 10 4 GeV 2 




Figure 15: Percentage difference in the NNLO up and down valence quark PDFs at Q 2 = 10 4 
GeV 2 relative to central (as(M§) = 0.118) set for fits with different values of as, with the 
percentage error bands for the central set also shown. 


the various PDFs as a function of x for Q 2 = 10 4 GeV 2 - a value of Q 2 relevant to data from 
the LHC. In almost every case the changes in the PDFs for the coupling varied in the range 
0.116 < as(Mz) < 0.120 are well within the PDF uncertainty bounds. 

As expected, the gluon distribution for x < 0.1 is larger for as(M §) = 0.116 and smaller for 
as(Mz) = 0.120: a change which preserves the product asg, which approximately determines 
the evolution of F 2 (x, Q 2 ) with Q 2 at low x. This is the dominant constraint on the gluon, and 
a smaller low x gluon leads to a larger high-x gluon (and vice versa ) due to the momentum 
sum rule. The u and d PDFs have the opposite trend as ag(Mf) changes. At small x values 
this is a marginal effect, due to the interplay of a variety of competing elements. At high x 
the decreasing quark distribution with increasing as is due to the quicker evolution of quarks 
to lower x. The insensitivity of the strange quark PDF to variations of as(M f) at low x is 
partly just due to the relative insensitivity of all low-x quarks, but is also partially explained 
by the comments in the previous section about the MMHT analysis pi] of dimuon production 
in neutrino interactions - where the changes in ag(Mf) are, to some extent, compensated by 
changes in the B(D —>■ /i) branching ratio parameter. 

In Fig. 16 we compare the changes in the gluon PDF for different fixed values of ag(M|) at 
a much lower value of Q 2 , namely Q 2 = 10 GeV 2 . Here the gluon PDF is much more sensitive 
to the value of as(M f), and the changes in the gluon PDF lie outside its uncertainty bounds. 
The message is clear. At the high value of Q 2 = 10 4 GeV 2 the long evolution length means 
that the gluon PDF in the relevant broad x range about x ~ 0.01 is determined by PDFs at 
larger x, and is relatively insensitive to the parameters of the starting distributions. 


7 Benchmark cross sections 

In this section we show uncertainties for cross sections at the Tevatron, and for 7 TeV and 
14 TeV at the LHC. Uncertainties for 8 TeV and 13 TeV will be very similar to those at 7 TeV 
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Figure 16: Percentage difference in the NNLO gluon PDFs at Q 2 = 10 GeV 2 relative to central 
(ag(M§) = 0.118) set for fits with different values of as, with the percentage error bands for 
the central set also shown. 


and 14 TeV respectively. We calculate the cross sections for W and Z boson, Higgs boson via 
gluon-gluon fusion and top-quark pair production. 

We calculate the PDF and as(M §) uncertainties for the MMHT2014 PDFs |4] at the de¬ 
fault values of a's(Mf). We use a value of Aag(Mf) = 0.001 as an example, simply because 
PDF sets are readily available with ag(M|) changes in units of 0.001. However, for values 
similar to Aag(M|) = 0.001 a linear scaling of the uncertainty can be applied to a very good 
approximation. As explained in Section [HJ the full PDF+ag(Mf) uncertainty may then be 
obtained by adding the two uncertainties in quadrature. 

To calculate the cross section we use the same procedure as was used in [!]. That is, for 
W, Z and Higgs production we use the code provided by W. J. Stirling, based on the calculation 
in [IS], [20] and [21], and for top pair production we use the procedure and code of [22]. Here 
our primary aim is not to present definitive predictions or to compare in detail to other PDF 
sets, as both these results are frequently provided in the literature with very specific choices of 
codes, scales and parameters which may differ from those used here. Rather, our main objective 
is to illustrate the procedure for estimating realistic PDF+cts(Mf) uncertainties. 

7.1 W and Z production 

We begin with the predictions for the W and Z production cross sections. The results at NNLO 
are shown in Table [l] In this case the cross sections contain zeroth-order contributions in as, 
with positive NLO corrections of about 20%, and much smaller NNLO contributions. Hence 
a smaller than 1% change in ag(Mf) will only directly increase the cross section by a small 
fraction of a percent. The PDF uncertainties on the cross sections are 2% at the Tevatron and 
slightly smaller at the LHC - the lower beam energy at the Tevatron meaning the cross sections 
have more contribution from higher x where PDF uncertainties increase. The as uncertainty 
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cr 

PDF unc. 

as unc. 

W Tevatron (1.96 TeV) 
Z Tevatron (1.96 TeV) 

2.782 

0.2559 

+0.056 f +2.0% \ 
-0.056 V— 2.0%/ 

+0.0052 /+2.0%\ 
-0.0046 \— 1.8%/ 

+0.018 /+0.65%\ 
-0.020 V— 0.72%/ 

+0.0015 / +0.59%\ 
-0.0018 V— 0.70%/ 

W + LHC (7 TeV) 

W~ LHC (7 TeV) 

Z LHC (7 TeV) 

6.197 

4.306 

0.9638 

+0.103 /+1.7%\ 
-0.092 V—1.5%/ 

+0.067 / +1.6% \ 
-0.076 v—1.8%/ 

+0.014 / +1.5% \ 
-0.013 V— 1.3%/ 

+0.058 /+0.94%\ 
-0.065 V —1.0% ) 

+0.043 /+1.0%\ 
-0.043 1-1.0%; 

+0.0091 / +0.94%5 

-0.010 V —1.0% ) 

W+ LHC (14 TeV) 

W~ LHC (14 TeV) 

Z LHC (14 TeV) 

12.48 

9.32 

2.065 

+0.22 /+1.8%\ 
-0.18 V —1.4%/ 

+0.15 / +1.6% 5 
-0.14 V —1.5%/ 

+0.035 / +1.7%\ 
-0.030 v—1.5%/ 

+0.12 /+0.97%\ 
-0.14 1-1.1% ) 

+0.098 /+1.1%\ 
-0.11 1-1.2%; 

+0.020 /+0.97%\ 
-0.025 V —1.2% ) 


Table 1: Predictions for W ± and Z cross sections (in nb), including leptonic branching, obtained 
with the NNLO MMHT2014 parton sets. The PDF and as uncertainties are also shown, where the 
as uncertainty corresponds to a variation of ±0.001 around its central value. The full PDF±a,s(Mf) 
uncertainty is obtained by adding these two uncertainties in quadrature, as explained in Section [5} 


is small, about 0.6% at the Tevatron and close to 1% at the LHC, being slightly larger at 
14 TeV than at 7 TeV. Hence, the as uncertainty is small, but more than the small fraction 
of a percent expected from the direct change in the cross section with as- In fact the main 
increase in cross sections with as is due to the change in the PDFs with the coupling, rather 
than its direct effect on the cross section. From Fig. 14 we see that the up and down quark 
PDFs increase with as below x ~ 0.1 — 0.2 due to increased speed of evolution. From Fig. 13 


we note that the strange quark PDF increases a little with as at all x values. As already 
mentioned the Tevatron cross sections are more sensitive to the liigh-x quarks, which decrease 
with increasing as, so this introduces a certain amount of anti-correlation of the cross section 
with as- However, the main contribution is from low enough x that the distributions increase 
with as, so the net effect is an increase with as a little larger than that coming directly from 
the a s dependence of the cross section. As the energy increases at the LHC the contributing 
quarks move on average to lower x and the increase of the cross section with as increases - 
very slightly more so at 14 TeV than at 7 TeV. However, even at 14 TeV the total PDF±as 
uncertainty obtained by adding the two contributions in quadrature, is only a maximum of 
about 25% greater (for W~) than the PDF uncertainty alone if A as(M^) = 0.001 is used. 


7.2 Top-quark pair production 

In Table [2] we show the analogous results for the top-quark pair production cross section. At 
the Tevatron the PDFs are probed in the region x ~ 0.4/1.96 ~ 0.2, and the main production is 
from the qq channel. As we saw, the quark distributions are reasonably insensitive to as(Mf) 
in this region of x, as it is the approximate pivot point of the PDFs. Hence, there is only a 
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cr 

PDF unc. 

as unc. 

tt Tevatron (1.96 TeV) 
tt LHC (7 TeV) 
tt LHC (14 TeV) 

7.51 

175.9 

969.9 

+0.21 (+2.8%\ 
-0.20 V— 2.7%/ 

+3.9 / +2.2%^ 
-5.5 V— 3.1%/ 

+ 16 / +1.6% \ 
-20 V— 2.1%/ 

+0.17 (+2.3%\ 
-0.15 V— 2.1%/ 

+4.1 /+2.3%\ 
-3.3 v—1.9%/ 

+ 16 / +1.6% \ 
-14 V —1.4%/ 


Tabic 2: Predictions for tt cross sections (in nb), obtained with the NNLO MMHT2014 parton sets. 
The PDF and as uncertainties are also shown, where the as uncertainty corresponds to a variation 
of ±0.001 around its central value. The full PDF±ais(Mf) uncertainty is obtained by adding these 
two uncertainties in quadrature, as explained in Section [5] 


small change in cross section dne to changes in the PDFs with as- However, the cross section 
for tt production begins at order a|, and there is a significant positive higher-order correction 
at NLO and still an appreciable one at NNLO. Therefore, a change in as a little lower than 
1% should give a direct change in the cross section of about 2%. This is roughly the change 
that is observed. This is compared to a PDF only uncertainty of nearly 3% due to sensitivity 
to higher x quarks that occurs for W, Z production. 


At the LHC the dominant production at higher energies (and with a proton-proton rather 
than proton-antiproton collider) is gluon-gluon fusion, with the central x value probed being 
x ~ 0.4/7 ~ 0.06 at 7 TeV, and x ~ 0.4/14 « 0.03 at 14 TeV. As seen from the left plot of 
Fig. 13 the gluon decreases with increasing ag(Mf) below x = 0.1 and the maximum decrease 
is for x ~ 0.02 — 0.03. The as(M §) uncertainty on a t j for 7 TeV is about 2%, almost as large 
as at the Tevatron, with the gluon above the pivot point still contributing considerably to the 
cross section, so the indirect a,j(M§) uncertainty due to PDF variation largely cancels. For 
14 TeV the lower x probed means that most contribution is below the pivot point and there 
is some anti-correlation between the direct as variation and the indirect, with a reduced as 
uncertainty of 1.5%. At this highest energy the PDF only uncertainty has also reduced to about 
2% due to the decreased sensitivity to the uncertainty in high-x PDFs, the gluon in this case. 
At the Tevatron and for 7 TeV at the LHC the as (Mg) uncertainty is a little smaller than the 
PDF uncertainty, and the total is about 1.3 times the PDF uncertainty alone. At 14 TeV they 
are very similar in size, so the total uncertainty, for Aas(Mf) = 0.001 is about y/2 that of the 
PDF uncertainty. 


7.3 Higgs boson production 

In Table [3] we show the uncertainties in the rate of Higgs boson production from gluon-gluon 
fusion. Again, the cross section starts at order a 2 s and there are large positive NLO and NNLO 
contributions. Hence, changes in as of about 1% would be expected to lead to direct changes in 
the cross section of about 3%. However, even at the Tevatron the dominant x range probed, i.e. 
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cr 

PDF unc. 

as unc. 

Higgs Tevatron (1.96 TeV) 
Higgs LHC (7 TeV) 

Higgs LHC (14 TeV) 

0.874 

14.56 

47.69 

+0.024 / +2.7%\ 
-0.030 V—3.4%/ 

+0.21 / +1.4% \ 
-0.29 V —2.0%/ 

+0.63 / +1.3%\ 

-0.88 V —1.8%/ 

+0.022 /+2.5%\ 
-0.018 V—2.1%/ 

+0.23 /+1.6%\ 
-0.22 V —1.5%/ 

+0.71 7+1.5%\ 
-0.70 V —1.5%/ 


Tabic 3: Predictions for the Higgs boson cross sections (in nb), obtained with the NNLO MMHT 
2014 parton sets. The PDF and as uncertainties are also shown, where the as uncertainty cor¬ 
responds to a variation of ±0.001 around its central value. The full PDF±cts(M§) uncertainty is 
obtained by adding these two uncertainties in quadrature, as explained in Section [5] 



0 

PDF unc. 

as unc. 

Higgs Tevatron (1.96 TeV) 
Higgs LHC (7 TeV) 

Higgs LHC (14 TeV) 

0.644 

11.28 

37.63 

+0.021 7+3.3%\ 
-0.022 V— 3.4%/ 

+0.21 7+1.9%\ 
-0.20 V —1.8%/ 

+0.67 7+1.8%\ 
-0.59 V-l.6%; 

+0.011 7+1.7%\ 

-0.0088 V-1.4%1 

+0.15 7+1.3%\ 
-0.14 V —1.2%/ 

+0.51 7+1.4%\ 
-0.50 V —1.3%/ 


Tabic 4: Predictions for Higgs Boson cross sections (in nb), obtained with the NLO MMHT 2014 
parton sets.The PDF and a s are shown, with the a s uncertainty corresponding to a variation of 
±0.001 around the central value (as(Mf) = 0.120). The full PDF±cts(Mf) uncertainty is obtained 
by adding these two uncertainties in quadrature, as explained in Section [5] 

x ~ 0.125/1.96 ~ 0.06, corresponds to a region where the gluon distribution falls with increasing 
as(Mz) and at the LHC where x ~ 0.01 — 0.02 at central rapidity the anti-correlation between 
as(Mg) and the gluon distribution is near its maximum. Hence, at the Tevatron the total 
as(Mz) uncertainty is a little less than the direct value at a little more than 2%, and at the 
LHC it is reduced to 1.5%. In the former case this is a little less than the PDF uncertainty 
of ~ 3%, with some sensitivity to the relatively poorly constrained high-x gluon, while at the 
LHC the PDF uncertainty is much reduced due to the smaller x probed, and is similar to the 
as(M f) uncertainty. Hence for Aas(M§) = 0.001 the uncertainty on the Higgs boson cross 
section from gluon-gluon fusion is about \/2 that of the PDF uncertainty alone. 

We also repeat the study at NLO for the Higgs cross section. The results are shown in 
Tables |4] and [5] for the central values of as(Mf) = 0.120 and as(Mf) = 0.118 respectively. The 
uncertainties are very different in the two cases, with the central values of the cross sections 
being about 3% lower for as(M§) = 0.118 than for as(Mf) = 0.120. Both sets of predictions 
are about 30% lower than at NNLO, highlighting the large NNLO correction for this process. 
The PDF uncertainties are very similar to those at NNLO, though a little larger in detail. 
However, the as(M§) uncertainties are noticeably reduced, as the large variation in the NNLO 
(0(d!g)) cross section with as is now absent. 
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a 

PDF unc. 

as unc. 

Higgs Tevatron (1.96 TeV) 
Higgs LHC (7 TeV) 

Higgs LHC (14 TeV) 

0.625 

11.00 

36.61 

+0.022 /+3.5%\ 
-0.019 V—3.0%/ 

+0.21 /+1.9%\ 
-0.19 V —1.7%/ 

+0.65 / +1.8%\ 
-0.56 V —1.5%/ 

+0.0096 /+1.5%\ 
-0.0086 V —1.4%/ 

+0.15 /+1.4%\ 
-0.14 V —1.3%/ 
+0.51 /+1.4% \ 
-0.50 V —1.4%/ 


Table 5: Predictions for Higgs Boson cross sections (in nb), obtained with the NLO MMHT 2014 
parton sets.The PDF and a s are shown, with the a s uncertainty corresponding to a variation of 
±0.001 around the central value (ag(Mf) = 0.118). The full PDF±as(Mf) uncertainty is obtained 
by adding these two uncertainties in quadrature, as explained in Section [5] 

8 Conclusions 

The PDFs determined from global fits to deep-inelastic and related hard-scattering data are 
highly correlated to the value of as(M^) used, and any changes in the values of a:g(M§) must 
be accompanied by changes in the PDFs such that the optimum fit to data is still obtained. 
In [1] we produced PDF and uncertainty eigenvector sets for specific values of as(M|), guided 
by the values obtained when it was left as a free parameter in the fit. In this article we explicitly 
present PDF sets and the global fit quality at NLO and NNLO for a wide variety of ag(Mf) 

values, i.e. cts(Mf) = 0.108 to ag(Mf) = 0.128 in steps of Aa;g(Mf) = 0.001. Hence, we 

illustrate in more detail the origin of our best fit ag(M|) values of 

NLO: crg(Mf) = 0.1201 ± 0.0015 (68% C.L.), (7) 

NNLO: «s(Mf) = 0.1172 ± 0.0013 (68% C.L.), (8) 

already presented in [4j, but also present the uncertainties. We show the variation of the fit 
quality with a?g(Mf) of each data set, within the context of the global fit, and see which are the 
more and less constraining sets, and which prefer higher and lower values. We see that most 
data sets show a systematic trend of preferring a slightly lower a'g(Mf) value at NNLO than at 
NLO, but note that no particular type of data strongly prefers a high or low value of ag(M|). 
HERA and Tevatron data tend to prefer higher values, but are not the most constraining data. 
There are examples of fixed target DIS data which prefer either high or low values and similarly 
for the LHC data sets, which are new compared to our previous analysis ra- Indeed our best 
values of ag(Mf) are almost unchanged from ag(Mf) = 0.1202 (NLO) and «s(M|) = 0.1171 
(NNLO). They are also very similar to the values obtained by NNPDF of ag(M|) = 0.1191 
(NLO)[23] and ag(M|) = 0.1173 (NNLO) [21]. However, our extraction disagrees with the 
recent value «g(M|) = 0.1132 (NNLO) in [25]. We find agreement at the level of one sigma 
or less with the world average value of ccg(M|) = 0.1187 ± 0.0005, and this improves when we 
include the world average (without the DIS determinations included) as a data point in our fit, 
when we obtain as(M §) = 0.1195 (NLO) and ccg(M|) = 0.1178 (NNLO). Hence, our NNLO 
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value including as(M^) as an external constraint is in excellent agreement with the preferred 
value, as(Mz) = 0.118, for which eigenvector sets are made available. The PDF sets obtained 
at the 21 different values of as(M^) at NLO and NNLO can be found at [26] and are available 
from the LHAPDF library |2?’]. They should be useful in studies of as(M^) by other groups. 

In order to calculate the PDF+«s(Mf ) uncertainty we now advocate the approach pioneered 
in [18] of treating PDFs with as(M§) ± Aas(M§) as an extra eigenvector set. As shown in [IS] , 
provided certain conditions are met (at least approximately), the as (Mg) uncertainty may be 
correctly added to the PDF uncertainty by simply adding in quadrature the variation of any 
quantity under a change in coupling A as(Mf) as long as the change in as{M f) is accompanied 
by the appropriate change in PDFs required by the global fit. As examples, we have calculated 
the total cross sections for the production of W, Z, top quark pairs and Higgs bosons at 
the Tevatron and LHC. For W and Z production, where the LO subprocess is 0(a Q s ) and is 
quark-initiated, the combined “PDF+cts” uncertainty is not much larger than the PDF-only 
uncertainty with a fixed as- However, the additional uncertainty due to as is more important 
for top quark pair production and Higgs boson production via gluon-gluon fusion, since the 
LO subprocess now is 0(a|), though the details depend on the correlation between as(M^) 
and the contributing PDFs. 

In addition, we note that for any particular process the details of the uncertainty can now 
be explicitly calculated in a straightforward way using the PDFs we have provided in this 
paper, together with the procedure for combining PDF and as(Mf) uncertainty discussed in 
Section [5] 

Moreover, it is also straightforward to apply the procedure to determine the uncertainties 
coming from combinations of PDF sets obtained by global analyses of different groups. Using 
techniques given in [28] [29? 30] [3T] it is possible to combine different PDF sets at a preferred 
value of as(Mz) such that the central value and the uncertainty of the combination is correctly 
obtained. The procedure to determine the uncertainty due to variations of as(M^) is as 
follows. If each group used in the combination also makes available sets of PDFs obtained 
by repeating their global fit^Jwith as(M|) ± Aas(Mf), then an additional pair of PDF sets 
representing the a$(M f) variation of the combination can be obtained just by taking the average 
of the PDFs from each group obtained at as(M|) + Aas(M^), and by taking the average at 
as(Mz) — Aas(Mf). As a result the PDF+as(M§) uncertainty for any quantity calculated 
using the combined set is just the PDF induced uncertainty at the preferred value of as(M \) 
added in quadrature to the as(Mf) uncertainty determined from the two combined sets defined 
at as(Mz) ± Act's (M§). Hence, a user may determine for any process the optimum prediction, 
the PDF uncertainty, the as(M§) uncertainty and the complete PDF+a^Mf) uncertainty 
arising from the combination of a whole collection of different PDFs. 

5 For instance, if ag(M|) = 0.118 is the preferred value then repeating global fits at as(M§) = 0.117 and 
as(M^) = 0.119 would be sufficient to quantify the uncertainty due variations of as- 
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